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ABSTRACT 

We study contributions to the cosmic microwave background (CMB) bispectrum from non-Gaussianity 
induced by secondary anisotropics during reionization. Large-scale structure in the reionized epoch both 
gravitational lenses CMB photons and produces Doppler shifts in their temperature from scattering off 
electrons in infall. The resulting correlation is potentially observable through the CMB bispectrum. 
The second-order Ostriker-Vishniac also couples to a variety of linear secondary effects to produce a 
bispectrum. For the currently favored flat cosmological model with a low matter content and small 
optical depth in the reionized epoch t ;$ 0.3, however, these bispectrum contributions are well below 
the detection threshold of MAP and at or below that of Planck, given their cosmic and noise variance 
limitations. At the upper end of this range, they can serve as an extra source of noise for measurements 
with Planck of either primordial nongaussianity or that induced by the correlation of gravitational lensing 
with the integrated Sachs- Wolfe and the thermal Sunyaev-Zel'dovich effects. We include a discussion of 
the general properties of the CMB bispectrum, its configuration dependence for the various effects, and 
its computation in the Limber approximation and beyond. 

Subject headings: cosmic microwave background — cosmology: theory — large scale structure of 
universe — gravitational lensing 



1. INTRODUCTION 

The increase in sensitivity of upcoming cosmic mi- 
crowave background (CMB) experiments, especially satel- 
lite missions, raise the possibility that higher order corre- 
lations in the CMB temperature fluctuations, beyond the 



two-p oint function, may be experimentally detected and 
studi ed in detail to look for doviatione from Gauesianity. 



daily MAP|j and Planck surveyor^, with their higher sen- 
sitivity and angular resolution can in principle be used to 
study the bispectrum. It is therefore important to investi- 
gate what signals are expected in the currently favored cos- 
mological models an d how they might b e separate d from 
each other. Rec ently ^pergel fc Goldberg (1999) and Gold 



berg fc Spergel| (1999) studied the CMB bispectrum under 



In the absence deviations in the initial conditions, they 
may arise from the imprint of the non-linear growth of 
structures on secondary anisotropics. 



the context of non-Gaussian contribution due second or 
der gravitational effects, notably the correlation induced 
by gravitational lensing of CMB photons and seco ndary 



Ear ily theoretical work on the three-point correlation 
f unct i on wao aimed at diDtinguiohing between yariouo the 



anisotrop ics from the integrated Sachs- Wolfe (ISW; Sachi 



ories of the origin of th e fluct uatio ns (e.g ., |Falk et al 
1993; |Luo fc Schramm] 1993; ^u| 1994; Pangui et a" 
1994). Within the context of inflationary models, espe- 
cially within slow-roll inflation, this topi c has been further 
addre ssed in several recent pap ers (e.g., Gangui fc Martir] 
1999; " ■ "" " 



:op 
fc Wolfq 1967) effect and t he thermal Sunyaev-Zel'dovich 
(SZ; ISunyaev fc Zel'do^dch| 1980) effect 



Here, we present two additional effects resulting through 
reionization. The flrst one is the correlation induced by 
lensing through the secondary Doppler effect. The second 
one arises from the coupling of the second order Os triker 



Wang fc Kamionkowaki 1999). — In typical models. 



Visli niac effect (OV; Ostriker fc Vishniac 1986a, b; Vish 



the e xpected non Gaussian contribution ie below the level 
that can be detected due to cosmic variance limitations. 

Recently, there has been much renewed interest in the 
three point function and its Fourier analogue the bispec- 
trum, both in anticipation of the high precision satellite 
data and from analy ses of the COB E DM R data (e.g.. 



Hinshaw et a] 1995; Ferreira et aL 1998; Pando et al 



199 8). The dete ction reported by Ferreira et al. (1998) 



and Pando et al. (1998) has been subsequently shown to be 
associat ed with a know n systematic e rror in the data time- 
stream ( Banday et al. 1999; see also Bromley fc Tegmark 
1999). This should be taken as a cautionary tale that even 
if an experiment has the raw sensitivity to detect expected 
bispectrum signals, removing systematic effects to the re- 
quired level will prove challenging indeed. 

Nonetheless, the upcoming CMB experiments, espe- 



niacj 1987) to other linear secondary effects such as the 
ISW, Doppler and SZ effects. The resulting bispectrum 
is proportional to the square of the matter density power 
spectra. We also study the conflguration dependence of 
these effects to address means by which they may be sepa- 
rated from each other and any primordial nongaussianity. 
The layout of the paper is as follows. In § g, we re- 
view the background material relevant for understanding 
the CMB bispectrum and its statistical properties in the 
context of the adiabatic cold dark matter (CDM) models. 
A summary of useful properties of the Wigner-3j sym- 
bol is presented in the Appendix. In § ||, we detail the 
coupling between gravitational lensing angular excursions 
and secondary anisotropics. Results for the currently fa- 
vored cosmology are presented in § 0. In § pi we treat the 

^ http://map.nasa. gsfc.gov 
2http://astro.estec.esa.nl/Planck/; also, ESA D/SCI(6)3. 



coupling between Ostriker-Vishniac (OV) effect and sec- 
ondary temperature fluctuations; results for the favored 
cosmological model are presented in §^. We conclude in 
§0- 

2. PRELIMINARIES 

We first review the properties of adiabatic CDM models 
relevant to the present calculations. We then discuss the 
general properties of the angular bispectrum of the CMB. 

2.1. Adiabatic CDM Model 

The expansion rate for adiabatic CDM cosmological 
models with a cosmological constant is 



H' 



Hi [r!„,(l + zf + nK{i + zf + riAJ , (1) 



where Hq can be written as the inverse Hubble distance 
today Hq^ = 2997.9ft.~^Mpc. We follow the conventions 
that in units of the critical density SHq/SttG, the contri- 
bution of each component is denoted f2i, i = c for the 
CDM, b for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities firn = ^c + ^b and 
fix = 1 — J2i ^i' which represent the matter density and 
the contribution of spatial curvature to the expansion rate 
respectively. 

Convenient measures of distance and time include the 
conformal distance (or lookback time) from the observer 
at redshift z = 



^{z) = 



dz' 



and the analogous angular diameter distance 



- H-io-i/2 



1/2^ 



dA = H-'nj^''smh{Hon\^r) 



(2) 



(3) 



Note that as ^Ik —<■ 0, dA —^ r and we define r{z — cxj) = 

ro- 

The adiabatic CDM model possesses a power spectrum 
of density fluctuations 

(SikySik')) = i2nfS^ik - k')^A^{k) , (4) 

where 6^ is the Dirac delta function^ and 



A^(fc) = '5^(^ 



n+3 



T\k)., 



(5) 



in linear perturbati on theory. We use the fitting formulae 
of Eiscnstcin fc Hu| (1999) in evaluating the transfer func- 
tion T{k) for CDM models. Here, 5h is the amplitude of 
present-day density fluctuations at the Hubble scale; we 
adopt the COBE normalization for 5h ( Bunn fc White| 
1997). 

The density field may be scaled backwards to higher 
redshift by the use o f the gr owth function G(z), where 
5{k,r) = G{r)5{k,{)) ( peeblc^ 1980) 



'We also use this symbol for Kronecker deltas. 



(6) 



Note that in the matter dominated epoch G on a 



(1 



it is therefore convenient to define an auxiliary 



quantity F = G/a. 

The cosmological Poisson equation relates the density 
field to the fluctuations in the gravitational potential 



.^a„(^ 



1 + ^§^K 



F{r)S{k,0). (7) 



Likewise, the continuity equation relates the density and 
velocity fields via. 



-iGS{k,0) — 



(8) 



where overdots represent derivatives with respect to radial 
distance r. 

For fluctuation spectra and growth rates of interest here, 
reionization of the universe is expected to occur rather 
late Zri ^ 50 such that the reionized media is optically 
thin to Thomson scattering of CMB photons r ^ 1 . The 
probability of last scattering within dr of r (the visibility 
function) is 



re-^ = XHoTHil + zfer^ 



(9) 



Here T{r) = L drr is the optical depth out to r, X is the 
ionization fraction, 



TH = 0M91{1 -Yp)nbh, 



(10) 



is the optical depth to Thomson scattering to the Hub- 
ble distance today, assuming full hydrogen ionization with 
primordial helium fraction of Yp. Note that the ionization 
fraction can exceed unity: X = {1 — 31^/4)/(l — Yp) for 
singly ionized helium, X — {1 — Yp/2)/{l — Yp) for fully 
ionized helium. We assume that 



X{z) = l- -ei-ic (^I^ 
2 \V2Az 



(11) 



such that hydrogen reionizes smoothly but promptly at Zri! 
some aspects of the Doppler effect, in particular its power 
spectrum and to a lesser extent Doppler-Doppler-OV cou- 
pling, are sensitive to the sharpness of this transition. 

Although we maintain generality in all derivations, we 
illustrate our results with the currently favored ACDM 
cosmological model. The parameters for this model are 
nc = 0.30, Ub = 0.05, Qa = 0.65, h = 0.65, Yp = 0.24, 
n = 1, and 6h = 4.2 x 10~^. This model has mass fluctua- 
tions on the 8h Mpc^^ scale in accord with the abundance 
of galaxy clusters as — 0.86. A reasonable value here is im- 
portant since the bispectrum is nonlinearly dependent on 
the amplitude of the density field. We consider reioniza- 
tion redshifts in the range 5 < Zji ^ 40 or 0.025 ^ r < 0.5 
and assume Az/(1 + z) — 0.1. 

2.2. Bispectrum 

The bispectrum Bi-^i^i^ is the spherical harmonic trans- 
form of the three-point correlation function just as the 
angular power spectrum Ci is the transform of the two- 
point function. In terms of the multipole moments of the 
temperature fiuctuation field r(n). 



air-, 



j dnT{h)Yr{n) 



(12) 



the two point correlation function is given by 

C(fi,m) = (T(n)r(m)) 



E 



K™,«;.m.)i^r*(A)V'(rfi)-(i3) 

Iimil2m2 

Under the assumption that the temperature field is statis- 
tically isotropic, the correlation is independent of m 



{(^Inm'^hmi) — "hl2^mim2^h i 



(14) 



and called the angular power spectrum. Likewise the three 
point correlation function is given by 

i3(n,m,i) = (T(fi)T(m)T(i)) (15) 

^ J2{ai,„nai2m2aHm,)Yi;^'{n)Y,^\rh)Y^^^(\) , 

where the sum is over (Zi, mi), (^2, m.2), (^3,7713). Statistical 
isotropy again allows us to express the correlation in terms 
an ?7i-independent function. 



\(^limi ^l2'm2^l^m,^) 



h h h 

TOi 7712 ^,3 



B, 



hhh 



(16) 



Here the quantity in parenthese s is the Wigner-3j' symbol 
Its orthonormality relation Eq. (A2) implies 



B 



hhi-s 



I h h h 



^-^ \ 77ll 77i2 7713 
mirn2rn3 ^ 



\^limi^l2m,2^l3m,3/ -V '/ 



The angular bispectrum, Bi-^i^i^, contains all the infor- 
mation available in the three-point correlation function. 
For ex ample, the skew ness, the collapsed three-point func- 
tion of Hinshaw et aj (1995) an d the equilateral configura- 
tion statistic of Fcrrcira ct al. (1998) can all be expr essed 



as linear combinations of the bispectrum terms (see Gan- 
gui et 



al 1994 for explicit expressions). 



It is also useful to note its relation to the bispectrum de- 
fined on a small flat section of the sky. In the flat sky ap- 
proximation, the spherical polar coordinates {9, (p) are re- 
placed with radial coordinates on a plane (r = 2 sin6'/2 « 
0, 0). The Fourier variable conjugate to these coordinates 
is a 2D vector 1 of length I and azimuthal angle (pi . The ex- 
pansion coefficients of the Fourier transform of a given 1 is 
a weighted su m over 777 of th e spherical harmonic moments 
of the same / ( |White et al] 1999) 



«(1) 



so that 



47r ^ 

2/ + 1 ^ ' 



27r 



^hn^ 



(18) 



(a*(li)a(l2))-^C^,a5], 



,im(<j}i^~(t)i^) 



Likewise the 2D bispectrum is defined as 

(a(li)a(l2)a(l3)) = (2^)'<5°(li + I2 + l3)S(li, I2, I3) 



(19) 
(20) 



{hhhY'^ 



B 



hhi-s 



E 



xe 



"I>l3) 



„jmi(0ij- 



h 



[ h h 

\mi 77l2 — 77li — 7^2 



The triangle inequality of the Wigner-3j symbol becomes a 
triangle equality relating the 2D vectors. The implication 
is that the triplet (^1,^2, ^3) can be considered to contribute 
to the triangle configuration li, 12,13 = — li -I- I2 where the 
multipole number is taken as the length of the vector. 

2.3. Fisher Matrix 

To quantify the amount of information contained about 
a parameter pi in the bispectrum, we shall define the 
Fisher matrix Fij. The Cramer- Rao inequality (Kendall & 
Stuart 1969) says that the variance of an unbiased estima- 
tor oi Pi cannot be less than (F"^)^^. In terms of the like- 
lihood L of observing bispectrum elements Ba = Bi-^i^i^ 
(arranged as a data vector) given the true parameters p 
(called the fiducial model), the Fisher matrix is defined as 



a2lnL(B;p) 
dpidpj 



(21) 



Under the approximation that the likelihood is Gaussian, 
this expression becomes 



^■'^^ dp, ^ '""' dp, 



(22) 



The covariance matrix between the a and a' bispectrum 
term is in general a complicated quantity to calculate (see, 
Heavensl 1998). However there are two simplifying assump- 
tions that make this expression tractable. Firstly, since the 
CMB is expected to be nearly Gaussian, the dominant con- 
tribution to the covariance comes from the 6-point func- 
tion of the Gaussian field. This can be expressed in terms 
the power spectrum of all contributions to the field com- 
bined, i.e. the cosmic signal, detector noise, and residual 
foregrounds in the map, 



Cr^ =.Ci+ Q"°'^° + CI 



foreg 



(23) 



Secondly, assuming all-sky coverage^ the power spectrum 
covaria nce i s diagonal in /, 771 and the covariance then be- 
comes ( Luo 1994) 



Caa' = {Bhl2l3Bl'J'^l'^) 

= Cf°'Cf°'Cf°' 16^(123) + <5°(231) + <5°(312) 
+6^(213) + 5°(132) + 5°(321)1 , (24) 



where 



S^{abc)^Sr,,jr2A 



(25) 



Here, we have assumed that h ^ h, h i^ h, h "^ h- The 
covariance increases by a factor of 2 for two Vs equal and 
a factor of 6 for all 3 Vs equal. 

Under these assumptions, the Fisher matrix reduces to 



l3>l2>h 



9Biii2i3 dBi^ 



I2I3 



dpi dpj 



(26) 



■^Even satellite missions will not have all-sky coverage due to the 
need to remove galactic contamination. This will increase the co- 
variance of the estimators but we neglect such subtleties here. 



where 



Table 1 



hhh 



i'^tot y^tot i^-ftot 

^1 ^2 ^3 



(27) 



-'(2(3 



€1. 



Ki.€is] 



Note that the covariance between permutations of (/i, I2, I3) 
restricts the sum to I3 > h ^ h- 

Since the signal is expected to be smaU, in this paper 
we will be mainly concerned with the overall observability 
of the bispectrum. Consider the most optimistic scenario 
where the form of the bispectrum is considered known 
and the only parameter of interest is its amplitude, i.e. 
B11I2I3 = ^-B/j^;2;3|fid where the true value oi A ~ 1. The 
Fisher matrix tells us that the variance of the measure- 
ments of A is no less than cr^(v4) = {F^^)aa or 



,2 _ 



a^iA) 



E 

l3>l2>ll 



B 



hhh 



hl2h 



(28) 



which corres ponds to the statistic introduced by ^perge^ 
& Go dberg (1999). We will often plot the contribution to 



X from a single I3, summed over li,l2 as dx /dl^ 

2.4. Detector Noise & Foregrounds 

Detector noise will degrade the sensitivity of a given ex- 
perime nt to t he bispectrum through Eq. (E3). As pointed 
out by Knox (1995), detector noise in a given frequency 
channel v can be treated as an additional sky signal with 
a power spectrum 



= (^) 



w 



-\vy 



e^(v)t(t+i) 



(29) 



if the experimental beam is Gaussian with width 9{v) in ra- 
dians (the full- width-half-maximum is given by FWHM ~ 
\S\n29iv)). Here the sensitivity measure w^'^iv) is de- 
fined as the noise variance per pixel times the pixel area 
in steradians. Modern experiments have many frequency 
channels with independent pixel noise. In this case we 
inverse variance weight the noise so that 



(30) 



y^noisc 



^ q-ise(^) 



where the sum is over the different frequency channels. 

We will use the MAP and Planck satellite specifications 
to illustrate the effects of detector noise (see Table ^. Note 
that this is somewhat optimistic since the multifrequency 
coverage of these experiments will have to be use d to re- 
move foreground contamination. Tegmark et al. (1999) 



have shown that given current expectations for foreground 
contributions, the increase in C*°' due to foreground re- 
moval is not expected to exceed '-^ 10%. Of course, the 
underlying premise that foregrounds add noise as a Gaus- 
sian random field is certainly incorrect at some level. It 
may well be that the bispectrum of the foregrounds is the 
limiting factor for observations of the cosmic bispectruni; 
a subject which lies beyond the scope of this paper. Fig. || 
displays the total power as a sum of the cosmic spectrum 
and inverse- variance weighted noise spectrum for the MAP 
and Planck satellites. 



CMB Experimental Specifications 



Experiment 


V 


FWHM 


lO'^AT/T 


MAP 


22 


56 


4.1 




30 


41 


5.7 




40 


28 


8.2 




60 


21 


11.0 




90 


13 


18.3 


Planck 


30 


33 


1.6 




44 


23 


2.4 




70 


14 


3.6 




100 


10 


4.3 




100 


10.7 


1.7 




143 


8.0 


2.0 




217 


5.5 


4.3 




353 


5.0 


14.4 




545 


5.0 


147 




857 


5.0 


6670 



NOTES. — Specifications used for MAP and Planek. Full width at 
half maxima (FWHM) of the beams are in arcminutes. iii~^" = 
AT X FWHM X 7r/10800. 



2.5. Secondary Anisotropics 

Secondary anisotropics in the CMB are those that are 
produced well after recombination at z ~ 1000. They 
generally fall into two broad classes: those arising from 
gravitational effects and those arising from Compton scat- 
tering. 

The gravitational effects are due to gravitational red- 
shift and lensing. The differential redshift effect from pho- 
tons climbing in and out of a time-varying gravitational 
potential along th e line of sight is called the integrated 
Sachs- Wolfe (ISW; [Sachs fc Wol^ 1967) effe ct while fluctu- 
ations are still linear, and the Rees-Sciama ( Rees fc Sciama 
1968) effect for the second-order and nonlinear contribu- 
tions. The ISW effect is important for low matter density 
universes $7™ < 1, where the gravitational potentials de- 
cay at low redshift, and contributes anisotropics on and 
above the scale of the horizon at the time of decay (see 
Fig. |l| calculated from Eq. Ip^). The Rees-Sciama effect is 



small in adiabatic CDM models (seljak 1996) and we will 
not consider it further here. 

Gravitational lensing of the photons by the intervening 
large-scale structure both redistributes power in multipole 
space and enhances it due power in the density pertur- 
bations. The most effective structures for lensing lie half 
way between the surface of recombination and the observer 
in comoving angular diameter distance. In the fiducial 
ACDM cosmology, this is at z ~ 3.3, but the growth of 
structure skews this to somewhat lower redshifts. In gen- 
eral, the efficiency of lensing is described by a broad bell 
shaped function between the source and the observer. The 
curve labeled "primary" in Fig. Q^ includes the gravita- 
tional lensing effect on th e power spectrum as calculated 



by the CMBFAST co de QSeljak 

Zaldarriaga fc Seljak 1997); lensing is responsi 



Zaldarriaga 1996 and 
Die for the 



power-law tail of anisotropics for I > 4000 here. 

Rescattering of the photons during the reionized epoch 
can both generate and erase anisotropics. The primary 
anisotropics are reduced as exp(— r) by scattering on small 



scales / <; 50. Density fluctuations can lead to optical 
depth variations and hence a patchy erasure of fluctua- 
tions. For the optically thin conditions considered here, 
this is a negligible effect. 

The bulk flow of the electrons that scatter the CMB 
photons leads to a Doppler effect. Its effect on the power 
spectrum peaks around the horizon at the scattering event 
projected on the sky today (see Fig. |l|). The contributions 
are strongly dependent on the optical depth in the reion- 
ized epoch; we take t — 0.1 as our flducial model but ex- 
plore optical depths up to r = 0.5. On scales smaller than 
the horizon at scattering, the contributions are mainly can- 
celed as photons scatter against the crests and troughs of 
the perturbation. As a result, the Doppler effect is moder- 
ately sensitive to how rapidly the universe reionizes since 
contributions from a sharp surface of reionization do not 
cancel. Following c urrent thinking t hat the universe reion- 
ized promptly (see Haiman fc Knox 1999 for a review), we 
take Az = 0.1(1 -I- zh) in Fig. K 

The Sunyaev-Zel'dovich (SZf punyaev fc ZePdovichl 1980) 
effect arises from the inverse-Compton scattering of CMB 
photons by hot electrons along the line of sight. This ef- 
fect has now been directly imag ed towards massive galaxy 



clusters (e.g., Carlstrom et al. 1996), where temperature 



of the scattering medium can reach as high as 10 keV pro 
ducing temperature changes in the CMB of order 1 mK. 
Here, we are interested in the SZ effect produced by large- 
scale structure in the general intergalactic medium (IGM) , 
where the exact calculation on its significance as a sec- 
ondary anisotropy is not directly possible due to the un- 
known distribution and clustering properties of baryonic 
gas and its temperature structure. Under certain simplify 
ing, albeit largely untested assumptions detailed in § pT5 



its power spectrum is given in Fig. H^ for both linear and 
nonlinear contributions. Given the untested assumptions, 
these spectra should be taken as provisional, even in their 
order of magnitude. Note that the SZ effect also bears a 
spectral signature that differs from the other effects. We 
have assumed observations in the Rayleigh- Jeans regime 
of the spectrum; an experiment such as Planck with sen- 
sitivity beyond the peak of the spectrum can separate out 
these contributions based on the spectral signature. 

The Ostriker-Vishniac effect arises from the second- 
order mo dulation of the Dopp ler eff ect by de nsity fluc- 
tuations ( pstriker fc Vishniac 1986a; Vishniac 1987). Its 
nonlinear anal ogue is the kinetic SZ effect from large-scale 
structure ( [Hrj 1999); we denote it OV'^"'' here to avoid 
confusion with the thermal SZ effect. Due to its density 
weighting, the OV effect peaks at small scales: arcminutes 
for ACDM. For a fully ionized universe, contributions are 
broadly distributed in redshift so that the power spectra 
are moderately dependent on the optical depth t, which 
we generally take to be 0.1. 

Variations in the ionization fraction of the gas around 
the ep och of reionizat ion ca n also modulate the D oppler 



effect (Aghanim et al. 1996; Gruzinov & Hu 1998; Knox 
et al.[] 1998). This effect depends on the detailed phys- 
ical processes occurring at reionization though it proba- 
bly contributes substantially to the sub-arcminute scale 
anisotropy. Because of the small scale of the effect and its 
unknown relationship to the density field and hence the 
other secondary effects, we do not consider its contribu- 
tion to the power spectrum or bispectrum here. 




10-° 
10" 

10-^° 
^10-" 

1 10-^' 

10-^' 
10-" 
10-^= 



Fig. 1. — Power spectrum for the temperature anisotropies in the 
fiducialjVCDM model with t = 0.1 {(zh> = 13 ) Az = 0.1(1 -I- zh) 
(see § p.5| for details). The curve labeled "primary" actually includes 
the small ISW, Doppler, and lensing contributions. Note that the 
predictions for the SZ power spectrum are highly uncertain and 
frequency dependent. We have also shown the instrumental noise 
contribution of MAP and Planck calculated using parameters in 
Table 1 which is important for signal-to-noise calculations. 



In summary, for ACDM cosmologies with optical depths 
T ~ 0.1, the ISW contributions dominate at low Z- values, 
the Doppler contributions at intermediate Z-values, and 
the SZ contributions at small scales. Given the crudeness 
of the approximation for the SZ effect, one cannot rule out 
the possibility that OV contributions dominate at high-L 
The power spectra of these effects imply that the intrinsic 
contributions of the ISW, Doppler and SZ effects are com- 
parable but appear at different angular scales and arise 
from different redshifts. It is for this reason that we shall 
consider the bispectrum contributions from these effects 
in the following sections. 

3. LENSING EFFECTS: DERIVATION 

The general derivation of CMB bispectrum due to the 
coupling of le nsing with secondary anisotropies, originally 
c ons idered by Goldberg & Spergel (1999), is reviewed in 
It is applied to Doppler secondary anisotropies in 
In § 3.3 and 3.4 for comparison, we briefly revisit the 



3T 
3.2 



coupling between gravitation al lensing and ISW and SZ 
effects respectively studied by Goldberg fc Spergel (1999). 



3.1. General Considerations 

Large-scale structure deflects CMB photons in transit 
from the last scattering surface. These structures also give 
rise to secondary anisotropies. The result is a correlation 
between the temperature fluctuations and deflection an- 
gles. This effect cannot be seen in the two point func- 
tion since gravitational lensing preserves surface bright- 
ness: deflections only alter the temperature fleld on the 
sky in the presence of intrinsic, primary, anisotropies in 
the unlensed distribution. The lowest order contribution 
thus comes from the three-point function or bispectrum. 

In weak gravitational lensing, the deflection angle on the 
sky is given by the angular gradient of the lensing potential 



which is itself a projection of the gravitational potential 
(see e.g. Kaiser 1992), 



e(m) = -2 r dr '!^^\^ /"\ ci>(r, mr) . (31) 

This quantity is simply related to the more familiar con- 
vergence 

1 



«(m) = -V^e(m) 



(32) 



^ , dA(i')dA(ro — r) ^^ ^, ^ , 
dr ^ '. ; -V]_<^{r, mr) , 



dA{ro) 



where note that the 2D Laplacian operating on $ is a spa- 
tial and not an angular Laplacian. The two terms k and O 
contain superficial differences in their radial and wavenum- 
ber weights which we shall see cancel in the appropriate 
Limber approximation. In particular, their spherical har- 
monic moments are simply proportional 



e. 



1(1 + 1) 



Ulr. 



dnl^"*(n)e(n) 



— <5(k)rr*(k)/r(fc) 



(33) 



with 



/r(fc)-/ drW'-''{k,r)Ukr), 



W'"'{k,r) = -mrn[^] F{r 



dA{rQ~r) 



(34) 



k J ^ "' dA{r)dA{ra) 
Here, we have used the Rayleigh expansion of a plane wave 
g.k. A. ^ 4^ ^ ,' J, (fcr)y ™* (k)y^" (fi) , (35) 

Ira 

and the fact that V^y," - -l{l + 1)^," — In an open 



the 



nnive^se^ one simply replaces the spherical Ressel functions 
with ultraspherical Bessel functions in expressions such as 
Eq. (g^. Since this is the case, we present the derivation 
in flat space. 

Since expressions of the type in Eq. ( p4[ ) frequently occur 
in the following calculations, it is useful to note that it 
may be evaluated efficientl y with a temporal version of 
the Limber approximatio n (Limber! 19 54) called the weak 
coupling approximation (Hu fc White 1996), 



lf{k)= drW''{k,r)j,{kr) 

Jo 



W 



W'^ik^l/k) drjiikr) (fc^ > 1) 

Jq W 



W^{k,l/k)^ 



A r[(? + i)/2] 

2k T[{1 + 2)/2] ' 



(36) 



and note that the ratio of gamma functions goes to y^jl 
for I 3> 1. Fo r the open universe generalization of this 
result, see Hu (1999). We employ this approximation for 
sufficiently high / values (usually I > 200), where the dif- 
ference between full integral and the approximation is suf- 
ficiently low (< 1%). 



Assuming a relation between the multipole moments 
and I^ as in Eq. (P3|), we find that the cross-correlation 
power between X and Y then becomes 



( 






^)^inJfA'{k)lHkW(k) 



(37) 



¥r^^^-^"(^''')^"(i'^)^^(^)' 



where the approximation involves a change in variables 
dA = l/k and we have restored generality for open ge- 
ometries where r ^ dA (see Hu 1999). Note that in the 
Limber approximation where dA is interchangeable with 
l/k, the superficial difference in the weights for k and 
disappears. 

The quantity of interest is the correlation between the 
deflection potential and secondary anisotropics 

TS(n)=^aLl'r(fi), (38) 

which becomes 

(e(n)rS(m)) = ^ (eLaL) Yr{n)Yr{rh) . (39) 

hn 

Again statistical isotropy guarantees that we may write 
the correlation as 



{^Im'^lm} =^1 = 



-2 






1(1 + 1) 

/dk 
-A2(fc)/f(fc)/;-(A), (40) 

«^/ drd^iyS(f,r)H^'-(f,r)A^(-l), 
l^ Jo dA dA dA 



where we have used eq uation (p3|) to relate the power 
spectrum bf defined by Goldberg & Spergel (199 9) and 
the K-secondary cross power spectrum defined by Seljak 



& Zaldarriaga (1999). The last line represents the Limber 



approximation and we have assumed that the secondary 
anisotropics are linearly related to the density field pro- 
jected along the hne of sight. 



d^k 



6{k)if{k)Yrik), 



/f(fc)= / drW^{k,r)ji{kr). 



(41) 



Note that the power spectrum of the secondary effect is 
then given by 



Cf = in 



dk 



A^{k)lf{k)lf{k). 



(42) 



We will discuss explicit forms for W^ for specific secondary 
effects in the following sections. 

Unfortunately, the secondary-lensi ng correlation is not 
directly observable. As pointed out by Goldberg fc Spergej 
(1999), it does however have an effect on the bispectrum 
which is in principle observable. The lensed temperature 
fluctuation in a given direction is the sum of the pri- 
mary fluctuation in a different direction plus the secondary 



anisotropy 



or 



~I][(aL+«L)V'(n)+ar. 

Im 

xVe(n)-Vy,'"(n) 

P S \ ^ p 

a/m = 0.1m + ^Im + 2^ '^/'rn' 

X /dnY,"*(n)Ve(n) • Vr,r'(n) 



(43) 



(44) 



Utilizing the definition of the bispectrum in Eq. (ITh, we 
obtain 



Bhhh 



E 



mini2m3 



h h h 
mi 1712 rns 



xVY(';"i*(m) • (Ve(m)T'5(n)) +Pcrm. (45) 

where the five permutations are with respect to the order- 
ing of {h.hjs)- 

Integrating by par ts and siniphfying further following 
Goldberg fc Spergel| (1999) leads to a bispectrum of the 
form: 



^l^2^3^ /(2;i + i)(2Z2 + i)(2;3 + i) 





An 
Gi6,3 + Perm. 



(46) 



where we have employed Eq. (A5) to perforin the angular 
integration. 

3.2. Doppler-Lensing Effect 

The Doppler effect generates temperature fluctuations 
as 

pro 

T'*°P(n) == / drg{r)h ■ v(r, fir) . (47) 

With the help of Eq. (|), 



47r 



h-k = J2^kYr{n)Yr{^). 



(48) 



and the Clebsh-Gordan coefficients for the addition of an- 
gular momenta with Z = 1, we can write the multipole 
moments as 

af°P = -4W I ^^y™*(k) 



(27r)3 



X / drgGj[{kr). 
Jo 



(49) 



Through integration by parts, this expression can be 
brought into the standard form of Eq. mW with 



W^°nk,r)^^{gG + gG), 



(50) 




Fig. 2. — The Doppler-lensing effect. Shown are the combined 
Doppler and double scattering effects, {solid line), the Doppler effect 
{dotted line), and the Limber approximation to the Doppler effect 
{dot- dashed line). At sufficiently high I, the difference between these 
three treatments can be ignored for most practical purposes. Top 
panel — The correlation power spectrum. Bottom panel — Contri- 
bution to x^ per log interval in (3. 



through integration by parts. 

Employing Eq. (Bfl) for the Icnsing correlation, we ob- 
tain 



, dop 



6tt' 



f 
xF(r) 



a„Ml I drd\{gG + gG) 



dAiro-r)^,^^^^^^ 



dAiro) 
in the Limber approximation for / > 200. 

3.3. Double- Scattering-Lensing Effect 



(51) 



As pointed out by Kaiser (1984), double scattering ef- 
fects must also be considered since the single scattering 
contributions are canceled along the line of sight at short 
wavelengths. Equation (49) implies the Doppler effect pro- 
duces an isotropic temperature fluctuation in the photons 
which subsequently last scatter into the line of sight. 



dop^ 

^00 



(fc)ro" = fc-^5'^°(fc,r)(5(k) 



J(k) 



dr'gGn[k{r~r')], (52) 



where g is to be interpreted as the visibility function for an 
observer at r, i.e. g = fexp[— T(r') -\-T{r)]. For scales that 
are much smaller than the width of the visibility function, 
i.e. approximately k{ro — r) 3> 1, one can take gG out of 
the integral leaving 



5'"°(fc,r) 



-G. 



(53) 



The contribution to the anisotropy today takes the stan- 
dard form of Eq. (|l]) 



In the Limber Umit, Eq. (Eo) can be siniphfied as 
67r2 



(54) 



br^ 



xF{r) 
again vaUd for £ ^ 200. 



ro 

VlmHl / drd\{gTG) 



dA(ro) 



(55) 



3.4. ISW-Lensing Ejfect 



The integrated Sachs- Wolfe effect ([Sachs fc Wolfe] 1967) 
results from the late time decay of gravitational potential 
fluctuations. The resulting temperature fluctuations in the 
CMB can be written as 



r^*^(n) 



dr^{r, hr) , 



(56) 



Using the Poisson equation (Eq. [Q), we can then bring 
the contributions into the standard form of Eq. ( |4l| ) 



Tyisw(fc) = _3f^^^^ ( E^ 



Hr)- 



(57) 



The two-point function produced between gravitational 
lensing angular deflections an d ISW effect can now be writ- 
ten as in the Limber limit as ([Goldberg fc Spergel 1999)^ 



bl 



ISW , 






^,nHo 



drFFd\ 



dAJrp -r) 
dAiro) 



A^l/dA). 



(58) 



For calculational purposes, we employ numerical integra- 
tion to / of 200 and use Limber approximated formulation 
thereafter. 

3.5. SZ- Lensing Effect 

The SZ effect leads to an effective temperature fluctua- 
tion in the Rayleigh- Jeans part of the CMB given by the 
integrated pressure fluctuation along the line of sight: 



T 



SZ/ 



drg 






5. 






(r, fir) . 



(59) 

where fc^ is the Boltzmann constant. Unfortunately, the 
clustering properties of the gas and its temperature struc- 
ture is uncertain even with state of the art hydrodynamic 
simulations. 

To obta in an order of magnit ude estimate of the effect, 
we follow poldberg fc Spergel (1999) in making several 
simplifying but largely untested assumptions. We assume 
that the gas is a biased tracer of the dark matter density 
(^gas = bgns^ and that the IGM temperature distribution 

^Note that our Limber approximated rnrrpla tinn hntwepn Ifrising 
and ISW effect is factor 2 lower than Eq. (14) of [Goldberg fc Spcrgej 
(1999); the same equation also contains an additional misprint witn 
respect to weighing by their r factors. 




Fig. 3. — Dependence of the Doppler-lensing effect on the ioniza- 
tion optical depth for t of 0.1, 0.3 and 0.5 with Az = 0.1(1 -I- ^ri)- 
Top panel — The 6; term for the Doppler-lensing effect. Middle 
panel — CMB power spectrum used in the noise calculation. Lower 
panel — The contribution to x^ P^r log interval in Z3 with MAP and 
Planck detector noise included. 



varies as T(,{r) — a(r)Tco. Finally we ignore temperature 
fluctuations in the gas. 

The effect can then be expressed as a weighted projec- 
tion of the density field of the form in Eq. (|4l|) 





W^'^ik) = -2A^^gaG , 




(60) 


ith a normalization given by 
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^ ksTco 










= 0.00783 fV^ 
V 4 J 


VlkeVy 




(61) 


or illustration 


3urposes, we ado 


Pt &gas = 


4 


and Too = 



IkeV throughout. 

The cross power between lensing and SZ can be written 
using Eq. (40) 



kSZ 



12^2 



2 aSZ 



-Hj^HqA 



drG'd\g 



dAiro 



dAiro) 



-A^ 



I 
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(62) 
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Fig. 4. — Linear and nonlinear density power spectra for the dark 
matter under the Peacock & Dodds (1996) scaling approximation for 
our fiducial ACDM cosmological model evaluated at the present. 



Since the SZ -lensing results in a correlation that peaks at 
/ ~ 100 (see, Goldberg & Spergel 1999), with no significant 
low Is 



contribution at low la, the Limber approximation may be 
used to calculate the whole effect. 

Both SZ effect and gravitational lensing angular deflec- 
tion potentials will be enhanced by nonlinear fluctuations 
in the density field. To get a qualitative understanding of 
the effects, assume that the gas density and potential fluc- 
tuations continue to track the dark matter in the nonlinear 
regime. The nonlinear evolution of the dark matter den- 
sity distribution has been well-studied with N-body simu- 
lations for adiabatic CDM cosmological models of interest 



here. We employ the scaling relation of Peacock fc Doddi 
(1996) to obtain the nonlinear power spectrum today (see 
Fig. 0) as well as its evolution in time. In Fig. n^, the curve 
labeled "SZ*^"'^" is the power spectrum formed by the SZ 
effect using the non-linear matter density power spectrum. 
The nonlinear effects generally increases the power due to 
SZ temperature fluctuations by a factor of 3 to 4 when I 
is in the interested range of 1000 to 5000. 

With the introduction of the non- linear power spectrum, 
we can write the cross-correlation between lensing and SZ 
as 



^SZ(nl) 



~1^ 



n, 



.HlA^^ 



drd\g 



dAiro) 



^2(nl) 



d^' 



(63) 



Since the nonlinear power spectrum no longer grows as 
G^, it must be evaluated along the line of sight at the 
corresponding lookback time. 

4. LENSING EFFECTS: RESULTS 

4.1. Doppler- Lensing Bispectrum 

The contribution to the bispectrum from the coupling of 
lensing to the Doppler effect is encapsulated in the power 
spectrum of the correlation b^ "^ (see Eq. [^ ) . Figure 
(top) shows its value for the Doppler effect in our fiducial 




Fig. 5. — Comparison of various lensing effects for our fiducial 
ACDM model and with t = 0.1 and Az = 0.1(1 -I- 2^). Top panel— 
The power spectrum of the correlation. Middle panel — Contribu- 
tions to x^ P6r log interval in Z3, assuming cosmic variance only. 
Lower panel — The same adding in detector noise for MAP and 
Planck. 



ACDM r = 0.1 cosmology assuming cosmic variance only 
C*°* = Ci (see Eq. ||). The zero crossing at I - 100 
is due to the competition between contributions on the 
reionization surface and those from the intermediate red- 
shifts. Likewise the Limber approximation to the integral 
also becomes excellent beyond I = 200. The inclusion of 
double scattering contributions to the Doppler effect only 
makes a minor difference for r ^ 0.1. 

We also show the contribution to the overall signal to 
noise squared (x^) per logarithmic interval in ^3 in Fig. 
(bottom). The minor differences in the x^ contributions 
at high ^3 between the Limber-approximated and Doppler- 
only spectra from the full calculation are primarily due to 
the differences of the 6; term at low I values. Low-/ values 
in bi contribute to high I3 values in x^ since the latter is 
a sum over h < h ^ h- The structure in the x^ plot and 
the rapid increase of its values around I ^ 1000 arises in 
large part from the structure of the primary C; itself which 
determines the Gaussian noise per mode. 

Since the signal is dominated by the smallest angular 
scales available in the measurements, it is important to in- 
clude the experimental beam and instrumental noise con- 
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tributions. In Fig. 0, we show the signal-to-noise ratio per 
mode of ^3 assuming MAP and Planck noise (see Table 
1). Here, we also show the affect of varying the optical 
depths from 0.1 to 0.5. Note that because reionization 
decreases the power spectrum of the anisotropics as e~^'^, 
the Gaussian noise in the bispectrum is also decreased with 
increasing optical depth. We have assumed a reionization 
width of Az = 0. 1(1-1- Zri) here, but the lensing correlation 
is not sensitive to this parameter. The lensing efficiency is 
a broad bell-shaped function in angular diameter distance 
and does not correlate well with contributions localized 
near the distant surface of reionization. 

As shown in Fig. 0, the bispectrum mode contribution 
to x^ for any Z3 due to coupling between reionized Doppler 
effect and gravitational lensing angular deflections is gen- 
erally less than 10~^, unless the optical depth to reion- 
ization is greater than what is currently inferred fro m the 



obser ved p ower spectr um of anisotropi cs (r < 0.5; Grif- 



fiths ^t al.| 1999; see, Haiman & Knox 1999 for a recent 
reviewj! Adding all the I3 mode contributions, the total 
X^ value is generally at a level below where one can expect 
its detection in the upcoming satellite data for reasonable 
optical depths (r ^ 0.3). 



SZ-lensing 



ISW lons i ng 



Doppler-lensing 



Fig. 



Configuration dependence of tiie secondary-lensing sig- 



nal. Plotted is the contributions to x of from single mode (3 = 300 
as a function of I2 and li . The configuration shapes are determined 
by bijth the power spectra of the cross-correlation (b; terms; see 
Fig. fi) and the CMB power spectrum. The spiky behavior in these 
plots IS an artifact of the plotting ipar the sharp crest introduced 
by the Wigner-3j symbol (see Fig. AlO). 



4.2. Comparison of Lensing Effects 

In Fig. g|, we compare the relative contribution to the 
CMB bispectrum from all the effects calculated. Both 
ISW and Doppler effects contribute at low I values to 
the correlation l^bi/2TT, while SZ effect to I ^ 100. The 
Doppler-lensing coupling is suppressed by a factor of ^ 8 
when compared to ISW-lcnsing coupling, which is in large 
part due to the inefficiency with which high redshift struc- 
tures lens the CMB. As shown, the SZ-lensing coupling 
and ISW-lensing coupling dominate the Doppler-lensing 
effect for optical depths t < 0.5. The nonlinear growth 
of fluctuations increases the SZ-lensing cross-correlation 
when I > 100 when compared to pure linear effects with 
a modest enhancement factor of ~ 5 at / '--^ 5000. Such 
an enhancement in the SZ-lensing cross-correlation power 
spectrum is consistent with increases in both SZ and lens- 
ing individual power spectra due to nonlinear effects. 

If untested assumptions on the SZ thermal effect were 
found to be false and the SZ contribution to be lower than 
currently presumed, than most of the contribution, at least 
in a low 0^ flat cosmological model, would come from 
coupling between gravitational lensing and the ISW ef- 
fect. However, even if this turns out to be the case, one 
may be able to use the spectral signature of the SZ ef- 
fect to separate its contribution. If SZ-lensing coupling is 
clearly detected, it would also support the current notion 
that most of baryons are present in the IGM as form of a 
gas at tem perature of ~ 1 keV. As suggested by Goldberg 
& Spergel (1999), t his would certain ly solve the "missing 
baryons" problem (Cen & Ostriker 1999). In the mean- 
time, improvements of our knowledge of the physical na- 
ture of the baryonic gas is clearly needed and can come 
from both observat ions as well as numerical simulations 
(e.g.,|daSilvaet alj 1999). 



Another obstacle for the unambiguous detection of the 
SZ, ISW and Doppler lensing effects is the ability to dis- 
tinguish their bispectrum contribution from t hose of other 



physica. 1 mechanisms. The x^ statistic of [Goldberg & 



Sperge] (1999) lumps all the information on each effect 



into one number and makes such an assessment difficult. 
Once all the effects are known the question of separabil- 
i ty c an be addressed with the Fisher matrix techniques of 
§2.3. In the meantime, it is useful to ask what combina- 
tion of (^1,^2,^3) is most of the signal coming from, i.e. 
what triangle configurations in the bispectrum dominate. 
In Fig. o, as an example, we show the general behavior of 
individual mode contribution to our signal-to-noise statis- 
tic x^ as a function of h and I2. Here, we have taken 
^3 = 300, which is representative of the their behavior in 
the ^3 range of interest. Since i?^(Zi, ^2, ^3) is zero when 
'1 + ^2 + ^3 is odd we have omitted these values for plotting 
purposes. 

The configuration dependences results from a combina- 
tion of the secondary-lensing cross-correlation power spec- 
tra and the CMB power spectrum. They are modulated 
by the behavior of the Wigner-3j symbol and variations 
in the noise term from the CMB power spectra As shown 
in Fig AlC, the Wigner-3j symbol also peaks at the lowest 
li that satisfies the triangle rule li — I3 — h- The CMB 
power spectrum noise introduces features associated with 
the acoustic peaks at I > 200 (see Fig. g and ||). 

The SZ-lensing bispectrum gets most of its contribution 
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when h ^ I3 — I2 with li taking values in the intermedi- 
ate I range of few hundreds where the cross-power peaks. 
This behavior also suggest why the SZ'-'^'^-lensing bispec- 
truni is not strongly enhanced even at high Z3 values by 
nonlinear effects, as most of the contributions come from 
intermediate 6; values, where nonlinear enhancement 

is negligible. 

For ISW-lensing, the cross-correlation term peaks to- 
wards low I values, and thus, its configuration dependence 
is such that the contributions to bispectrum peaks when 
/i ~ Z3 — I2 with ^1 in the range of few tens. This behavior 
continues out to high ^3 values. 

The Doppler-lensing configuration also shows a similar 
behavior but with contributions from somewhat smaller 
^1 ^ ^3 — h with Zi again taking values close to the peak 
in the cross-correlation power spectrum. These subtle dif- 
ferences in configuration may assist in isolating and iden- 
tifying the various contributions to the bispectrum once it 
is observed. 

5. OSTRIKER-VISHNIAC COUPLING: DERIVATION 

We consider the generation of CMB bispectra through 
coupling between the second order reionized Ostriker- 
Vishniac (OV) effect and two other linear sources of 
anisotropics. Since the OV effect is itself a secondary 
anisotropy source, we shall see that it only couples to 
other secondary anisotropy sources. First, we present a 
general derivation of the coupling between OV effect and 
secondary effects and in § 



5.1 



and then turn to specific sec- 
ondary effects such as the ISW, SZ and reionized Doppler 
effect. 

5.1. General Considerations 

The OV effect arises from the modulation of the Doppler 
effect by dens ity fluctuations whic h affect the prob ability 
of scattering (Ostriker & Vishniac 1986a; Vishniac 1987). 
The OV temperature fluctuations can be written as 

r (n) — I drg{r)n ■ v(r, nr)5{r, iir) 

d^k f d^k' 



= ^ / drgGG , ^^^^3 j ^^^^3 

X(5(k - k')5(k')e*-"'' [n • D(k, k')] , (64) 



where, 



D(k,k') 



1 



k-k' 

|k-k'|2 



k^ 



(65) 



is the projection vector that couples the density and ve- 
locity perturbations. 

We can now expand out the temperature perturbation, 
T*-*^, into spherical harmonics: 



„ov* _ ■ 



''Im 



dn / dr (gGG) 



d^ki 



xe 



(27r)3 J (27r)3 



r(ki)<5*(k2) 
(66) 



where we have symmetrizised by using ki and k2 to rep- 
resent k and k — k' respectively. Now the dot product 
between n and the projection vector is 



n-D(ki,k2) = 



1 /n • ki n • k2 



kf 



^2 



(67) 



with 



A-TT 

n.k = ^— fcrr(AW*(k) 



(68) 



Since the dot product is symmetric, we only consider one 
term here, and multiply the result by 2. By using the 
Rayleigh expansion (Eq. [p5|), we can rewrite the multi- 
pole moments as 
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.(4^)3 



dr 



d^ki f d^k; 
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limi l2'm2 fn' 
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xi^7^*(ki)>^r'*(ki)i^r*(k2) 

dnrr(n)y,7Hn)y,7(n)rr'(A)- 



(69) 



The power spectrum of the OV effect may be calculated 
from this expression. Since the end e xpre ssion is cumber- 
some we revert to the expressions of Hu (1999) in calcu- 
lating the OV power spectrum shown in Fig. n^. 

To leading order in A'^{k), bispectrum contributions in- 
volve one OV term and two linear sources of anisotropics. 
Recall that a general linear effect can be ex pre ssed as a 
weighted projection of the density held (see |3.l| ) 



''Im 



fk 

2^ 



6{k)lf{k)Yr{k), 
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(70) 



After some straightforward but tedious algebra, we can 
write 



„s s OV* 



(47r 



A2(fci)A2(fc2)/f,(fci) 



dki f dk2 

Xllik2)[I^l{k„k2)+I^l{k2,k,)] 

x|dnl^™^(n)y™^*(n)y™^*(n), (71) 



where 



I^l{ki,k2)= I rfrW^OVji,(fc2r)j^fcir), 
1 



VF"^(fci,r) 



ki 



gGG. 



(72) 



In simplifying the integrals involving spherical harmonics, 
we have made use of the properties of Clebsch-Gordon 
coefficients, in particular, those involving 1 = 1. 
In order to construct the bispectrum, note that 

(«'m,4m.«°™3) = i-^yHalraAm.a?^*^,) ■ (73) 

Under the assumption that "S" denotes the sum of all the 
sources so that the two contributions are indistinguishable, 
the bispectrum becomes 
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(74) 



where we have used Eq. ([A5|), Eq. (|A2|), and Eq. (|A6|). 
Here, 

y<I?ZAkiM)ll{ki)lUk2) ■ (75) 



Note that we have rewritten the ki -^ fc2 term in Eq. (71) 
as an li -^ I2 interchange so that in Eq. wm "Perm/' 
means a sum over the remaining 5 permutations of {li,l2,h) 
as usual. In the following sections, we evaluate this expres- 
sion with 



This Limber approximation reduces the dimension of the 
integrals from 5 to 1. Note that where the Limber approx- 
imation applies only equal time correlations contribute. 
On small angular scales then only secondary, not primary, 
anisotropics couple to the OV effect in the bispectrum. 

5.2. ISW-ISW-OV coupling 



Recall that from § 3.4, the weight function for the ISW 
effect is 



M^^=^(/fc,r) 
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(80) 



Substituting this weight function into Eq. ([79f) leads to 



.isw-isw_ 367r^ 2 tt4 r\rA2 
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(81) 



If{k) = //SW(^) ^ ^SZ(^) ^ ^dop(^) ^ ^SW(fc)^ (7g) 

The last term is the Sachs- Wolfe effect which we explicitly 
consider to show that the OV effect mainly couples with 
secondary not prima ry a nisotropies. We drop the double 
scattering effect of §3_^ as it is a small contribution for 
small optical depths. Note that expanding I^ in Eq. (75) 
produces many cross-terms which we will call hybrid ef- 
fects. 

In general, Eq. (|75|) involves five integrations, three over 
radial distances and two over wavenumbers. As in the lens- 
ing bispectrum calculation, these integrals can be simpli- 
fied using the Limber approximation for sufficiently large 
{li, ^2)- Here, we employ a version based on the complete- 
ness relation of spherical Bessel functions 



dkk'F{k)ji(kr)]i{kr') « -d-^S^^ir - r')F{k%_^ , 

2 dA 

where the assumption is that F{k) is a slowly- varying func- 
tion. Applying this to the integral over fc2 yields 

bl:iX^-rlJdrJdnJ^l,A^ik,)A\k2) 
W^{k2,r)W'^^{ki,r) 



xW^iki,n)- 

XJ/'i(fcl?')jii(fcl''l), 



r/2 



(78) 



Integrating by parts and assuming negligible boundary 
terms yields. 
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(79) 



We employ this Limber approximation in our calculation 
of the ISW-ISW-OV bispectrum effect. 

5.3. SZ-SZ-OY 

In a similar manner, we can ca lcul ate the SZ-SZ-OV 
effect following the discussion in j |3.5| where the weight 
function is given as 



VF^^(fc,r) 



-IJ^'^agG, 



(82) 



and is independent of the wave vector k. Now the Limber 
approximated bispectrum term, used in the calculations 
presented here, is 
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(83) 



Similar to our calculation on the enhancement of SZ- 
lensing bispectrum due to nonlinear growth of density fluc- 
tuations, we also consider the effect of nonlinearities on the 
SZ-SZ-OV bispectrum. As shown in Fig. H^, non-linearities 
enhance both SZ and OV effects and at high /, the OV 
effect involves the large-scale velocity field and small scale 
density field in the nonlinear regime (Hu 1999). For the 
bispectrum coupling, therefore, we replace only one of the 
power spectrum terms with the nonlinear relation, 
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(84) 



As before, with SZ-lcnsing nonlinear cross-correlation, 
the SZ-SZ'^"'^-OV coupling is calculated with the nonlinear 
density power spectrum evaluated along the line of sight 
at the corresponding lookback time. 

5.4. Doppler-Doppler-OV 
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The reionized Doppler effect was discussed in detail in 
T2| and has a weight function 



W'^°'P{k,r) 



1 

fc2 



(gG + gG). 



(85) 



Unhke ISW-ISW-OV and SZ-SZ-OV, where the Limber 
approximation is sufficiently accurate for our purposes, 
explicit integration over wavenumber is necessary here. 
The Limber approximation breaks down in the large-angle 
regime and hence affects the coupling between the veloc- 
ity field in the OV effect and the Doppler effect. On the 
other hand, the density field in the OV effect is dominated 
by small-scale fluctuations where the Limber approxima- 
tion is excellent. This implies that we may use Eq. ( |7q ) in 
calculating the Doppler-Doppler-OV effect: 
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dk 
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The integrals above include a derivative of the Bessel func- 
tion and is numerically difficult to evaluate. As in Eq. ([79|), 
one can integrate by parts to obtain a more tractable form. 
We considered both these approaches and the resulting bis- 
pectra agree at a level of 10%, with most of the difference 
resulting from the numerical computation of the coupling 
term involving the derivative of the Bessel function. This 
agreement also suggests that the evaluation at fc = l/dA in 
the Limber approximation can be performed after all the 
integration by parts are complete as is assumed in Eq. ([79|) . 

5.5. ISW-SZ-OV 

In addition to coupling between OV effect and secondary 
anisotropics of similar kind, we also consider hybrid cou- 
plings. The hybrids that contribute the most involve a 
large scale effect such as the ISW effect to couple with the 
velocity field in the OV effect, and a small-scale effect such 
as the SZ effect to couple with the density field in the OV 
effect. 
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To estimate the effects of nonlinearities, we take the ap- 
proximations introduced for the SZ-SZ("'^-OV above and 
find 
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5.6. Doppler-SZ-OV 



Similar to ISW-SZ-OV effect, the hybrid coupling of the 
reionized Doppler and SZ effects to OV effect. As with 
the Doppler-Doppler-OV calculation, we integrate the cou- 
pling between Doppler effect and OV velocity part and 



use the Limber approximation to describe the coupling 
between SZ and OV density part. The bispectrum term is 
then, 



y.G'g'aGdAfi,{kr)lf°^{k). 



(89) 



As with the other OV-SZ couplings, the nonlinear gener- 
alization of the effect is given by 



,dop-sz(ni) _ 167r- s^ f dk 
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xg'aGdAj'i,{kr)lf°^{k). (90) 

5.7. SW-SZ-OV 

Finally, we consider a hybrid involving the Sachs- Wolfe 
(SW) effect at the last scattering surface and the SZ effect. 
The weight function for the SW effect can be written as 

W^^ik,r.) = -^^FS^{r-r.), (91) 

and is evaluated at the last scattering surface, r» « rp. 
The Limber approximation implies that there is no cou- 
pling between the SW effect and the SZ or OV effects 
because there is no overlap in their weight functions. To 
the extent, that such a coupling exists the Limber approx- 
imation fails. 

Since it is at large scales that the Limber approximation 
breaks down, we are interested in the coupling between OV 
velocity part and the SW effect. As in the Doppler effects, 
we use Eq. ( [78|) to evaluate the deviations from the Limber 
approximation. 
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and its nonlinear analogue 
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Here, we do not consider couplings such as SW-SW-OV 
effect as the OV density part, with a redshift window at 
low redshifts, does not couple to SW effect at the last 
scattering surface. 

6. OSTRIKER-VISHNIAC COUPLINGS; RESULTS 

In Fig. ra, we show x^ contributions per logarithmic in- 
terval of ^3 for the coupling between linear secondary ef- 
fects and the OV effect as derived in the previous sec- 
tion. Here, we have assumed cosmic variance in the noise 
only. For the coupling between two secondaries of the same 
type, the bispectrum contributions per /3 are substantially 
smaller than their lensing counterparts (see Fig. ph. An 
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Fig. 7. — Contribution to x^ per log interval in Z3 for the OV coupling effects for the fiducial ACDM model with r = 0.1, A^ = 0.1(1 + Zri) 
and no instrumental noise. Le/i "panel — Coupling with two secondary effects of the same kind. Right panel — Hybrid effects involving 
secondaries of different kind. Other than Doppler-Dopplcr-OV and Doppler-SZ-OV couplings, the secondary-secondary-OV bispectra are only 
mildly sensitive to the reionized optical depth and width. 



SZ-SZ-OV 




Fig. 8. — Configuration dependence of the SZ-SZ-OV (top) and 
ISW-SZ-OV {bottom) couplings with (3 = 100. This figure is anal- 
ogous to the one shown in Fig. a for secondary-lensing coupling. 
Contributions generally peak towards (3 ~ I2 and low h for all ^3 
values. 



examination of the configuration dependence of the con- 
tributions shows that most of the contributions come from 
triplets which makes triangles with h <$^ h < h (see Fig. g 
for an example). The reason for this behavior is that the 
OV effect involves the large scale bulk velocity field and 
small scale density fiuctuations. One secondary effect cou- 
ples to the former and one to the latter creating the de- 
sired configuration. Because secondary effects tend to be 
strongly peaked to either small or large angular scales, any 
coupling involving two secondaries of the same type suffers 
suppression on one end or the other. 



Such configuration dependence also suggests that hybrid 
couplings should contribute more strongly. In particular 
in the last section, we considered hybrids composed of one 
large-scale secondary effect coupling to the OV velocity 
field and one small-scale secondary effect coupling to the 
OV density field. As shown in Fig. 0, the contribution 
to the bispectrum produced by the hybrid ISW-SZ cou- 
pling is substantially greater than either the ISW-ISW or 
SZ-SZ couplings. The same is true for the Doppler cou- 
plings: Doppler-SZ coupling yields a larger effect than ei- 
ther Doppler-Doppler or SZ-SZ coupling. It is however a 
smaller contribution than the ISW-SZ coupling. 

Except for couplings involving the Doppler effect, the 
results are fairly insensitive to the reionized optical depth 
since the others are weighted toward low redshifts. In- 
creasing the optical depth to 0.3 from 0.1, as shown in 
Fig. 7, the Doppler-Doppler-OV effect is increased by fac- 
tor of ~ 10, and when the optical depth is further increased 
to ~ 0.5, the Dopplcr-Dopplcr-OV contribution becomes 
comparable to that of the ISW-ISW-OV bispectrum sig- 
nal at ^3 of a few thousand. The hybrid Doppler-SZ-OV 
bispectrum is less sensitive to the reionized optical depth 
than the symmetric Doppler-Doppler-OV coupling, and is 
below the bispectrum produced by the hybrid ISW-SZ-OV 
for optical depths to reionization of current interest. 

Nonlinearities in the density field can enhance couplings 
involving the SZ effect. Under the simplifying assumption 
that the gas density traces the dark matter density, the 
enhancement in x^ for these effects is within a factor of 
100, at ^3 of few thousands, which is consistent with the 
enhancement in the SZ and OV power spectra due to non- 
linearities (see Fig. |^), and the behavior of li, I2 and ^3 con- 
figurations in producing the OV bispectrum. The smaller 
enhancement factor in the SZ'"')-lensing bispectrum, when 
compared to secondary-SZ'^'^-OV effects, comes from the 
fact that the secondary-lensing cross-correlation is not en- 
hanced by non-linear effects out to Z of '^ 100 (see. Fig. ||), 
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Fig. 9. — Contributions to the x'^ as a function of ^3 for our fiducial 
ACDM model for ISW-SZ-OV and its non-linear analogue ISW- 
SZ '"'-OV with MAP and Planck noise included in the variance. As 
shown, the signal-to-noisc ratios are generally below the detection 
level for MAP and near or below the threshold for Planck (x^ ~ 0.3,3 
for the linear and non-linear effects). These effects are the largest 
of the secondary-OV couplings considered. 



while the cross-correlation between SZ and density part 
of the OV effect is even enhanced at very low I values, as 
can be seen from Fig. ||. One should note that the current 
calculation using the non-linear power spectrum should be 
taken as an estimate of the upper limit of the effects since 
on the smallest scales gas pressure will make its distribu- 
ti on sm oother than that of the dark matter (see discussion 



Hu 1999). 



In general, Limber approximation when applied to cou- 
pling between OV effect and linear sources of anisotropics 
suggests that coupling only exists for effects with strong 
temporal overlap in the weight functions W^ and W'^^ . 
This suggests that OV coupling with all primary effects, 
which contribute at z ~ 1000 are heavily suppressed. We 
tested this by numerically integrating the coupling to to 
the SW effect. As shown in Fig. |^, the mismatch in red- 
shift between the OV effect and the SW leads to a large 
suppression of bispectrum signal and the contribution is 
generally the lowest of all secondary-secondary-OV cou- 
plings considered here. 

In Fig. g, we study the possible detection of OV- 
secondary couplings by future satellite CMB missions. 
Here, we have only focussed on the hybrid coupling of ISW 
and SZ effects to OV effect (ISW-SZ-OV), and its non- 
linear analogue ISW-SZ^^'^-OV, as they have the largest 
bispectrum signal of all secondary anisotropics that couple 
to OV effect. The signal-to-noise ratio generated by such 
effects is, however, small and their detection is marginal 
at best even with Planck. The ISW-SZ("')-OV effect has 
the largest x^ with a total value of ^ 2.8 for our fiducial 
ACDM cosmological model. In a flat universe with a cos- 
mological constant, for reasonable variations in Clm from 
0.2 to 0.7, more or less consistent with current observa- 
tional constraints on this parameter, the total x^ does not 
vary significantly to make this effect certainly detectable 
with Planck. The variation in x^ with ri,„ is such that the 



largest x^ value is ~ 3.1 when flm ^ 0.3 and falls below 
1 when rtm > 0.7 and tends to zero with Qm —^ 1; the 
fall-off is primarily due to the decrease in ISW effect with 
increasing $!,„ and the ISW effect is zero in an Einstein-de 
Sitter universe with ri,„ = 1. 

Beyond MAP and Planck, the largest contribution to 
the signal-to-noise for a perfect experiment comes from 
/3's corresponding to arcminute scales and are at a level 
that eventual detection may be possible. Even with such a 
high resolution experiment, one should bear in mind that 
their detection requires configurations in I space with one 
short side which picks up contributions from the ISW ef- 
fect. Its detection perhaps might involve a combination 
of a small angular scale experiment such as the planned 
ground-based interferometers and the satellite missions 
which have the larger angular coverage. The tight lo- 
calization of the bispectrum contributions of modulated 
Dopplcr effects such as the OV effect is also a useful prop- 
erty from the perspective of their role as a contamination 
for the secondary-lensing measurements or intrinsic non- 
gaussianity limits. Even if the amplitude is enhanced, say 
by contributions from the analogous patchy-reionization 
effect (see §^^) , a large part of the contamination may be 
removed by eliminating only a small range in the bispec- 
trunr terms. 

7. DISCUSSION 

Let us summarize the results of this study. Gravita- 
tional lensing angular excursions couple with first order 
secondary anisotropics generated at reionization to pro- 
duce a bispectrum in the CMB. Here, we have consider the 
coupling between lensing and the Doppler effect and have 
shown this to be a significantly smaller effect on the bis- 
pectrum than that between lensing and SZ and ISW effects 
for the low optical depths (r <, 0.3) expected in adiabatic 
CDM models. For the currently favored ACDM model, 
under the assumptions made here on the physical state 
of baryonic gas distribution, the bispectrum produced by 
SZ-gravitational lensing angular deflections is generally 



higher t han that produced by ISW- lensing (Goldberg 
Sperge] 1999). Results from higher resolution hydrody- 
namic simulations and constraints from observations are 
clearly needed to improve the calculation associated with 
SZ effect. 

The OV effect couples with two secondary sources of 
anisotropics to produce a bispectrum. Given that OV ef- 
fect is due to a product of density and velocity flelds, the 
bispectrum contributions are maximized when one of the 
secondaries peaks at large angular scales and the other at 
small angular scales. Such hybrid couplings produce rel- 
atively large signals in particular the ISW-SZ-OV effect 
if the assumptions involved in estimating the SZ contribu- 
tion prove to be correct. For reasonable models of reioniza- 
tion and currently favored cosmological models, however, 
the contribution to bispectrum by secondary-secondary- 
OV coupling is below the level that can be expected to 
be detected by MAP and is marginal at best for Planck. 
Although these signals are unlikely then to be detected 
in Planck, they do serve as a source of residual system- 
atic errors for other measurements of the bispectrum that 
can grow to be comparable to the cosmic variance term 
around the arcminute scale. Removal of these contribu- 
tions is facilitated by the strong localization of these ef- 
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fects at I3 ~ I2, li ^ 10. Removal of a broader range of li 
is undesireable since many of the interesting bispectrum 
effects are also maximized for I3 ^ h- 

In general, bispectrum contributions from the coupling 
of secondary anisotropics depends not only on the intrinsic 
amplitude of the secondary effects but also their overlap 
in redshift. In the small-scale Limber approximation, only 
equal time correlations contribute. Many of the effects 
considered here are substantially smaller than one would 
naively guess due to mismatch in the epochs at which the 
effects contribute their signal. 

We have considered contributions to the CMB bispec- 
trum from reionization to leading order in the density fluc- 
tuations. This is appropriate for the angular scales probed 
by the upcoming satellite missions. On arcminute scales 
and below, CMB anisotropics will be dominated by contri- 
butions from truly non-linear structures in the gas density, 
temperature and ionization state. The nongaussianity in- 
duced on these quantities by structure formation will be a 
rich field for future studies as experiments begin to probe 
the subarcminute regime of CMB anisotropics. 
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APPENDIX 

USEFUL PROPERTIES OF THE WIGNER-3J SYMBOL 

Here, we review the properties of the Wigner-3j symbol 
that are useful for the derivations in the text. It is defined 
by its relation to the Clebsch-Gordan coefhcient. 
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and as a consequence obeys the orthonormality relation 
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Fig. AlO. — Absolute value of the Wigner-3j symbol for m^ = as 
a function of ^i and I2 when I3 = 300. We plot only even Ii+l2 + l3- 



From the Clebsch-Gordan relation and series one can 
derive the integral relation to spherical harmonics. 
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From conjugation of this relation, one finds that 
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For TTii 
uated 



0, the Wigner-3j symbol can be efficiently eval- 
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for even X; it vanishes for odd L. In Fig. AlO, we show 



the absolute value of the Wigner 3-j symbol when mi — 
as given above. This expression is encountered in all 
our calculations on the bispectra produced by secondary 
effects. Note that other Wigner 3j terms, especially when 
rui is not equal to zero may be efficiently evaluated through 
recursion relations a nd the WKB approx imation at high 
multipolc argument ( Schulten fc Gordon] 1975). 
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where L = h + h + h- Given symmetry under spatial 
inversions only bispectrum terms with even L are non- 
zero. Since parity forces L to be even, all permutations 



and ni cgation of the ttl^s arc equal. Furthermore, the an - 
gular I momontum coloction ruloc require li < \ lj ^fc | for 



all permutations of the indices, and ttii + 7772 + "^3 — 0. 
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